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Computer simulations of a two-dimensional lattice 
of magnetic dipoles are performed on the Connection 
Machine. The lattice is a discrete model for thin dims 
of amorphous rare earth-transition metal alloys, 
which have application as the storage media in eras- 
able optical data storage systems. In these simula- 
tions the dipoles follow the dynamic Landau Lifsh- 
itz-Gilbert equation under the influence of an effec- 
tive held arising from local anisotropy, near- 
neighbor exchange, classical dipole-dipole interac- 
tions. and an externally applied held. Various sou- 
rces of coercivity, such as defects and or in- 
homogeneities in the lattice, are introduced and the 
subsequent motion of domain walls in response to 
external heids is investigated. 

Key words: coercivity. magnetization reversal 

l. Introduction 

Magnetization reversal in thin dims of amorphous 
rare earth-transition metal (RE-TM) alloys is of con- 
siderable importance in erasable optical data stor- 
age. 1 ' 131 The success of thermomagnetic recording 
and erasure depends on the reliable and repeatable 
reversal of magnetization in micron-size areas within 
the storage medium. A major factor usually en- 
countered in descriptions of the thermomagnetic 
process is the coercivity of the magnetic medium. 
Coercivity. however, is an til-defined concept general- 
ly used in the phenomenology of bulk reversal, but 
its usefulness and relevance to the phenomena occ- 
urring on the spatial/temporal scales of thermo- 
magnetic recording are seriously in doubt. To begin 
with, there is the problem of distinguishing nuclea- 
tion coercivity from the coercivity of wall motion. 
Then there is the question of speed and uniformity of 
motion as the wall expands beyond the site of its 
origination. Finally one must address the issues of 

J. Mag. Soc. Japan. Vol. 15. No. 1. 1991 


stability and erasability — also intimately related to 
coercivity —in a framework w ide enough to allow the 
consideration of local instability and partial erasure 
It is fair to say that the existing models and or 
theories of coercivity 14 31 are generally incapable of 
handling the problems associated with thermo- 
magnetic recording and erasure. In our view the 
natural vehicle for conducting theoretical investi- 
gations in this area is dynamic computer simulations 
based on the fundamental concepts and equations of 
micromagnetics. We have conducted several such 
studies in the past 32 ' -35 ' and have succeeded in clari- 
fying some of the issues involved in the nucleation 
stage of the magnetization reversal process. Also 
in.-. zed in these studies were the dynamics of wall 
motion in homogeneous media. The purpose of the 
present paper is to investigate the coercivity of wall 
motion at the submicrometer scale and in the pres- 
ence of defects and or inhomogeneities, using large 
scale computer simulations. 

The simulations were performed for a two dimen- 
sional hexagonal lattice of magnetic dipoles follow- 
ing the Landau-Lifshitz-Gilbert equation. 3 * 1 In addi- 
tion to interacting with an externally applied fieid. 
the dipoles were subject to effective fields arising 
from local uniaxial anisotropy, nearest neighbor ex- 
change. and long range dipole-dipole interactions. 
Details of the micromagnetic model have been previ- 
ously published 32 37 ~’ 39f and will not be repeated here. 
Suffice it to say that the massive parallelism of the 
Connection Machine on which these simulations 
were performed, together with the fast Fourier trans- 
form algorithm 37 381 which was used to compute the 
demagnetizing field, enabled us to accurately simu- 
late a large (256 - 256) hexagonal lattice of dipoles. 
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Since the lattice constant was chosen to be 10A in 
the^e simulations, the total area of the lattice corre- 
"ponits a section or the magnetic him with dirnen- 
t "Jon am. L y ^ JJJ a:n Periodic boundary 
o>nd::i"-:!* uere imposed ■ *n - he lattice m noth V and 
y directions. Out the film thickness along Z was 
finite and ’he magneti/Mtu «n distribution in this di- 
rection was . i spumed ' *■».■ '.in dorm. 

In a previous paper where nucleation ot reverse- 
magnetized domains was under consideration, it was 
>hnwn that the fields required to initiate the reversal 
process m a uniform material are generally higher 
than those observed in practice Various submicron- 
size detects” were then introduced in the magnetic 
state of the lattice and the values of nucleation coer- 
civity corresponding to different types, sizes, and 
strengths of these "defects” were computed, Voids, 
for instance, were found to have an insignificant 
effect on the value of the nucleation field, but re- 
verse-magnetized seeds, formed and stabilized in 
areas with large local anisotropy, could substan- 
tially reduce the nucleation coercivity. Similarly, 
the presence of spatial variations in the magnetic 
parameters of the material, such as random axis an- 
isotropy. was found to affect the coercivity of nucle- 
ation. 

Random spatial fluctuations and structural mag- 
netic defects also create barriers to domain wall 
motion. These barriers are overcome only when 
sufficiently large magnetic fields (in excess of the 
so-called wall motion coercivity) are applied. Simula- 
tions reveal that wall coercivity in amorphous 
RE-TM alloy films is generally less than the corre- 
sponding nucleation coercivity. This finding is in 
agreement with the experimentally observed square 
shape of the hysteresis loops in these media. The 
strength of wail coercivity. of course, depends on the 
tvpe and size of fluctuations and/or defects. The 
results reported in the next section are intended to 
clarify some of these relationships. 

Magnenzation-pattem-displays in this paper uti- 
lize a color coding scheme. Since the magnitude of 
the magnetization vector M will be fixed throughout 
the lattice, the Color Sphere is used to represent its 
local orientation. The Color Sphere is white at its 
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Fig. 1 The Color Circle is used to encode the 
direction of magnetization m the plane of the 
lattice. In this scheme a red pixel is associated 
with local magnetization direction along -.V. 
light green corresponds to - Y. blue to - X. and 
purple to - Y When a vector is not completely 
in the plane of the lattice, but has a perpendicu- 
lar component along -Z (or -Z). its associated 
color is obtained by mixing the color of its in- 
plane component with a certain amount of white 
(or black), the strength of white (or black- de- 
pending on the magnitude of the vertical com- 
ponent. A vector fully aligned with the -Z 
direction is shown by a white pixel, while a 
vector in the - Z direction is displayed as black. 

north pole, black at its south pole, and covers the 
visible spectrum on its equator in the manner shown 
in Fig. 1. As one moves from the equator to the north 
pole on a great circle, the color pales. i.e.. it mixes 
with increasing amounts of white, until it becomes 
white at the pole. Moving towards the south pole has 
the opposite effect as the color mixes with increasing 
amounts of black. Thus, when the magnetization 
vector at a given site is perpendicular to the plane of 
the lattice and along the positive (negative) Z axis, its 
corresponding pixel will be white (black). For \f in 
the plane of the lattice the pixel is red when pointing 
along ~X. light green along * Y, blue along - A. and 
purple along - Y. In the same manner, other orienta- 
tions of M map onto the corresponding color on the 
Color Sphere. 
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Fi* 2 Formation of domain walls in a medium with random axis anisotropy (cone angle 0-45*1 and in 
the absence of an applied field, (a) Dipoles in the white region are initialized along ~Z, while dipoles in the 
dark region are initialized along -Z. (b) The state of the lattice at t =0.8 ns. Each wall contains three 
vertical Bloch lines (VBLs) at this stage, (c) The state of the lattice at t- 1.7 ns. The number of V BLs has 
not changed since the previous frame, but they have moved along the walls, (d) The steady state of t e 
lattice at f = 4.56 ns. The number of VBLs in each wall is still three. 


2. Structure and Motion of Domain Walls 

The following set of parameters (used in all simula- 
tions reported in this paper) are typical of amorphous 
films of TbFeCo used in magneto-optical recording: 

saturation magnetization 100 emu/cm 2 3 * * * , aniso- 

tropy energy constant K u = 10* erg/cm 3 , exchange 

stiffness coefficient A t = 10' 7 erg/cm, film thickness 
h =500 A. damping coefficient a = 0.5, and gyromag- 
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netic ratio t— — 10' Hz/Oe. In subsequent discus- 
sions the term “random axis anisotropy is meant to 
imply that the anisotropy axes of the lattice are 
distributed randomly and independently among the 
lattice cells (or among various groups of these cells). 
By keeping the deviation angle 6 of the anisotropy 
axis from Z below a certain maximum angle & the 
random assignment of axes preserves the perpendic- 
ular nature of the overall anisotropy. Thus 0 is 
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Fig. 3 Collapse of the stripe domain depicted in Fig. 12td j under the external field H t 1000 Oe. Frames 
a; and «bi show the state of the lattice at f = 0.96 ns and f^3.58 ns. respectively. 



Fig. 4 Patchy lattice with 346 patches. Colors 
are used solely for the purpose of identification 
of different patches, otherwise they have no 
meaning or significance. 

selected with equal probability from the values in 
the interval [0, &]: for brevity. 0 will be referred to as 
the cone angle in this paper. No constraints are im- 
posed on the distribution of the azimuthal angle <p 
and its values are therefore selected (also randomly 
and independently) from the interval [0. 360 5 ]. It is 
perhaps worth mentioning at this point that the 
random assignment of the axes of anisotropy to indi- 
vidual dipoles does not automatically result in rapid 
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variations of the direction of magnetization in space. 
In fact, the strong exchange field in our simulations 
gives rise to a smooth distribution of the magnetiza- 
tion vector across the lattice, even w hen large cone 
angles are involved. 

Figure 2 shows the structure of domain walls in a 
uniform medium with the aforementioned set of pa- 
rameters. The lattice for this set of simulations had 
random axis anisotropy over individual cells with a 
cone angle of 0=45°. that is, the axis of local an- 
isotropy was selected randomly for each cell of the 
lattice with 0<=[ 0. 45’ ]. 0^[O. 360° ]. and the values of 
6?‘and 0 chosen for one cell were independent of those 
selected for all the other cells. Initially the central 
band of the lattice was magnetized in the -Z direc 
tion while the remaining part was magnetized in the 
-Z direction, as shown in Fig. 2(a). When the lattice 
was allowed to relax for 0.8 ns. the pattern in Fig 2 
(b) was obtained. Notice that there are three vertical 
Bloch lines (2?r VBLs) in each wall and that the walls 
are no longer straight. By allowing the lattice to 
relax for another 0,9 ns we obtained the pattern of 
Fig. 2(c). which shows significant VBL movements 
along the walls. Finally. Fig. 2(d) shows the steadv 
state situation at / = 4.56ns. Both walls are now 
straightened considerably, but the number of VBLs 
in each wall has not changed; no amount of relaxa- 
tion can unwind a 2,t Bloch line. 
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A perpendicular held // f 200 Oe. applied to the 

state of Fig. 2(d\ moves the walls somewhat closer to 
each other, but fails to eliminate them The force of 
demagnetization apparently opposes the external 
field in collapsing the reverse magnetized stripe. The 
-Jumain will collapse, however, under an applied field 
, 1000 Oe as shown in Fig 3 Frames 'a 1 and 
foj m thi" figure correspond, respectively, to t 0 96 
ns and t 3.58 ns. where t ts the elapsed time since 
the application of the field Although not shown in 
this figure, the walls proceeded to collide and an- 
nihilated each other at a later time. 

In the above example the randomness in the lattice 
parameters is clearly too weak to cause significant 
coercivity tor the wall motion. To remedy this prob- 
lem. w e increased the fluctuations* correlation length 
by introducing large patches, that is. patches whose 
dimensions are large compared to the width of the 
domain wail. Figure 4 shows a typical lattice cov- 
ered with 346 patches of random shape and size. 
These patches were created by selecting at random a 
number of lattice sites as seeds, and growing out- 
ward from them (in a random fashion) until every 
site in the lattice belonged to one patch or another. 
By assigning different attributes to different patches 
one can thus create spatial variations in the structure 
magnetic properties of the lattice over length ales 
comparable to the average patch size. 

Figure 5(a) shows a stripe of reverse magnetization 
in the patchy lattice of Fig. 4. For the sake of clarity, 
the boundaries of the patches are highlighted in the 
figure. Each patch is assigned an axis of anisotropy, 
randomly and independently of all the other patches, 
with a cone angle of 0=45 J The walls in Fig. 5 are 
somewhat more jagged than those in Fig. 2, where 
the patches were essentially the size of an individual 
lattice cell. Under an applied field of H, = - l 5 kOe 
the walls in Fig. 5(a) moved slightly and then came to 
equilibrium as shown in Fig. 5(b). A plot of the 

average magnetization of the lattice along theZ axis. 
versus time is shown in Fig. 5(c). This plot indicates 
that the initial <Af,X which is slightly above zero, 
has increased to about 20% of its saturation value m 
the first 0 5 ns after the application of the field, but 
has stopped growing at the point. Compare this 
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situation with that depicted m Figs. 2 and 3. w here • 
kOe of field was sufficient to annihilate the walls 
Clearly it ts the presence of patches and not the 
demagnetizing force: which is responsible for block 
ing the wall motion Thus the wall motion coercivity 
has increased as a result of increased correlation 
among the local easy axes. When the field was 
further increased to //< 2 kOe. it became possible to 

push the walls beyond the barriers and force 'hem to 
collide and disappear 1 Because ot the periodic noun 
dary condition, the walls collide after wrapping 
around the lattice.) Figure 5'd- shows an advanced 
state of domain expansion under the 2 kOe held t 
1.154 ns). Note that the lower wall has remained 
almost intact, while the upper wall has moved suo- 
stantially The plot of .V/, • versus time in Fig 5ie 
starting wuth the application of the 2 kOe fieid 
reveals that the movement is slow in the beginning, 
as the field struggles to overcome the pinning of the 
walls. Once released, the wall moves rapidly tor a 
period of about l ns until either another pinning 
occurs or the demagnetization force begins to push 
the walls apart (remember the periodic boundary 
condition). In any event, the motion slows down at 
/=l.2ns and the growth rate of .W, drops by 
nearly a factor of 2. The walls, however, continue to 
move out. wrap around at the boundary, and eventu- 
ally annihilate. 

In order to understand the effect of patch size on 
coercivity. we repeated the preceding simulation for 
another lattice which had the same set of parameters 
as the lattice in Fig. 5. but whose total number of 
patches had increased to 1300. Again we found that 
H t = l.SkOe could not move the walls significantly, 
but that H,-2 kOe could. It is probably safe, there- 
fore. to assert that the average size of the patch does 
not affect the coercivity in a substantial way. so long 
as it is larger than the characteristic width of the 
domain wall in the given medium. 

Figure 6(a) shows another strip of reverse mag- 
netization in the patchy lattice of Fig. 4. This time, 
however, a few of the patches have been made void 
by assigning the value of zero to their magnetic 
parameters. These void patches are shown as gray 
regions in the figure. (No special property is assumed 
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Fig. 5 Strip of reverse magnetization in the 
patchy lattice of Fig. 4. (a) The state of the 

lattice at //, = 0. (bl At t= 1.5 ns after the applica- 
tion of a 1.5 kOe field, ic) Average magnetization 
(M t ) versus time under the applied field of H t - 
1.5 kOe. (d) The state of the lattice 1.154 ns after 
the application of H t = 2 kOe. (el Average mag- 
netization (M t > versus time under the applied 
field of H t = 2 kOe. 
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Fig. 6 Strip of reverse magnetization in a 
patchy lattice with voids, (a) At H t — Q, (b) After 
t- 5 ns under the applied field of //, = 15 kOe. (c) 
Total magnetic energy versus time under H t = 1.5 
kOe. 


Fig. 7 Lattice with 662 patches, of which 24 
are “defective" A normal patch is assigned an- 
isotropy axis within a cone of 0=10° from the 
perpendicular direction, whereas the easy axis of 
a defective patch is nearly in the plane of the 
lattice. The color coding scheme depicts the 
direction of anisotropy for the individual patches. 
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Fig. 8 Magnetization distribution in the patchy lattice of Fig. 7. (a) The remanent state at H, = 0. (b) The 
steady state under H t - -8900 Oe. The remaining frames show the time evolution of the lattice under H t 
- -8920Oe. (c) At / = 0.24 ns the reversal is underway with a nucleus forming in the lower left corner. 
The periodic extension of this nucleus in the upper right corner of the lattice is also visible, <d) At f = 0.72 
ns. (e)Atf=lns. (f) At / = 1.08 ns. (g) At t= 1.248 ns. (h) At t = 1.32 ns. 


for the dipoles at the void boundaries, their magnetic 
parameters being the same as those elsewhere in the 
lattice. Of course no exchange field is exerted on the 
boundary dipoles from the neighboring cells on the 
void side, and the magnetic charges that accumulate 
on the void boundaries are automatically accounted 
for when the demagnetizing field is computed.) The 
remaining patches are all identical in their magnetic 
properties except for the value of the anisotropy 
constant which fluctuates randomly and inde- 
pendently from patch to patch. (The standard devia- 
tion of these fluctuations is 20% of the average value 
of K a .) No other spatial variations in the parameters 
have been assumed and, in particular, all axes of 
anisotropy are perpendicular ii.e.. 0=0). 

The walls in Fig. 6<a) have automatically adjusted 
themselves to minimize their lengths by attaching to 
the voids in the neigh borhood. Minimization of 

length is tantamount to minimization of total wall 
energy and is therefore favored by the magnetic 
system under consideration. Figure 6(b) shows the 
state of the lattice under an applied field of H t = 1 5 

kOe at t - 5 ns. Apparently, the walls have continued 
to seek voids to attach onto, while expanding in 


response to the external field. Figure 6(0 is a plot of 
the total magnetic energy of the lattice during this 
growth period. It is marked by slow declines, charac- 
teristic of continuous wall motion, and rapid drops, 
corresponding to detachments from or attachments 
to the voids. It is thus observed that void-like defects 
in real media can create jagged domain boundaries, 
increase the coercivity of wall motion, and cause 
discontinuities in the propagation process. 

As another example consider the case of a patchy 
lattice with a different kind of “defect”. The 662 
patches in this case were divided into two groups as 
shown in Fig. 7. The first group, consisting of 638 
patches, were assigned nearly-perpendicular axes of 
anisotropy, as has been our practice so far. Each one 
of the patches in this group was assigned an axis, 
randomly and independently of the others, with a 
cone angle of 0= 10 J The second group of 24 pat- 
ches were assigned axes of anisotropy that were 
nearly parallel to the plane of the lattice. The ran- 
domly and independently chosen values of Q and <p 
for these patches belonged to the intervals [80 J . 90 ] 
and [0. 360°]. respectively. (In the color-coded distri- 
bution of the anisotropy axes shown in Fig. 7, the 
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Fig. 9 E volution of the magnetization state of Fig. 8(e) after the removal of the external field, (a) f = 0.44 
ns. ib) t = 1.12 ns. (c) t = 1.16 ns. (d) f = 5.83 ns. 


patches with prominent colors are those with nearly 
in-plane axes.) The remaining parameters of the 
lattice were the same as before. In the following 
discussion this lattice shall be referred to as “the 
lattice with in-plane defects'*. 

The entire lattice with in-plane defects was initial- 
ly saturated along the i - Z axis and then allowed to 
relax in zero field, until the remanent state of mag- 
netization, shown in Fig. 8(a), was arrived at. The 
strong exchange interaction has clearly forced the 
magnetization of the in-plane patches towards the 
~Z axis, but the tendency towards the plane is still 
visible in some areas. We subjected the remanent 


state to a reverse-magnetizing external field along 
and searched for the critical magnitude of this 
field that would initiate the reversal process. The 
critical field for nucleation was found to be 8910 Oe 
It is important to note that in the absence of the 
in-plane patches the same sample had a nucleation 
coercivity close to 18kOe. The in-plane regions, 
therefore, facilitate the nucleation process. In Fig. 8 
(bt we show the pattern of magnetization distribu- 
tion under an applied field H t = -8900 0e. which is 
only slightly weaker than the critical field. The 
defects are seen to have been pushed toward the 
plane, yet the field is not quite strong enough to 
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cause them to reverse. When the field was raised to 
H t ~ - 8920 Oe. the reversal began, as the time evolu- 
tion series in Figs. 8(c) to 8(h) shows. First the defect 
at the lower left corner of the lattice (with its periodic 
extension at the upper right) nucleated. Then the 
nucleus grew until the entire sample was reverse- 
magnetized. Note how the defects seem to attract the 
domain wall as it approaches them, and then try to 
pin the wall to prevent its further progress. The 
applied field, however, is strong enough to overcome 
the pinning of the wall and bring the reversal to 
completion. A hysteresis loop with high squareness 
is the hallmark of this type of reversal. 

In order to investigate the pinning process and the 
phenomenon of wall coercivity, we took the state of 
the lattice in Fig. 8(e) as the initial state for another 
simulation and set the external field to zero. Note 
that because of the periodic boundary condition im- 
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Fig. 10 fa) Total magnetic energy £ T and aver- 
age magnetization (Xf t > of the lattice during the 
relaxation process shown in Fig 9. The total 
energy is the sum of exchange, anisotropy, and 
demagnetizing energies which are shown sepa- 
rately in (b) and !c). 


posed on the hexagonal lattice, the several black 
regions in Fig. 8(e) are in fact different pieces of one 
and the same reverse-magnetized domain. This 
(roughly circular) domain must begin to collapse 
immediately after the removal of the external field. 
The reason for the collapse is that the domain wall 
energy density of the material \o m =4 V AJ( U = 1.265 
erg cm 2 ) is greater than the demagnetizing energy 
density at the wall ( ^ — 6zr/iA/f =0.943 erg cm 2 ) and. 
therefore, energy minimization favors the shrinking 
of the domain. Figure 9 shows the time evolution of 
the shrinking process, with frames ( a ) to id) corre- 
sponding to / = 0.44 ns. 112 ns. 1.16 ns. and 5.83 ns 
after the external field is reset to zero. In the absence 
of defects the collapse would have been complete, 
but in the present situation after an initial period of 
shrinking, the defects trap the domain wall and pre- 
vent its collapse. Note in Fig. 9(c), for instance, how 
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Fig. 11 ia> Magnetization distribution pattern corresponding to the relaxation of the state in Fig. 8<fi 
under zero applied field. The steady state shown here is obtained after 4 ns of relaxation. <b) Dot product 
of the initial state of magnetization shown in Fig. 8(f) and the final (relaxed) state shown in Fig 1 ii a j. The 
color coding scheme for this figure which is different from that used elsewhere in the paper is described 
in the text. 


the defect in the lower central part of the lattice 
attracts the wall and keeps it pinned there after- 
wards. Therefore, like the other types of defect 
described in previous examples, defects with in-plane 
anisotropy are likely responsible for the observed 
wall coercivity of the RE-TM media. A comparison 
of Fig. 9<d) with Fig. 7 clearly indicates that the 
domain is fully anchored on a number of defects. 
The fact that the domain in Fig. 9(d) is stable has 
practical significance in data storage applications, 
since it provides one possible mechanism for the 
stability of very small thermomagnetically recorded 
marks. 

Figure 10 shows plots of energy £ and average 
magnetization \f t > during the relaxation process 
depicted in Fig. 9 As expected, the total energy 
of the system decreases in time, while the average 
magnetization along Z increases. A sudden jump in 
the curve of energy (such as the one at t- 1 14 ns) is 
an indicator of the capture of the wall by a defect. 
Plots of the various components of energy in Figs. 10 
(b) and 10(c) show that the exchange and anisotropy 
energies of the system generally decrease with time, 
while the demagnetizing energy increases. This is to 
be expected of course in the present situation since. 


as was pointed out earlier, during relaxation the total 
length of the domain wall must decrease. 

Finally, starting from the state of Fig 8(f). we show 
in Fig. 1 1 the result of relaxation under zero applied 
held. The magnetization pattern in Fig. 1 1(a) is the 
steady state of the lattice obtained after 4 ns of relax- 
ation. The reverse domain here is somewhat larger 
than the one in Fig. 9(d), simply because it started 
from a larger initial domain. The jaggedness of the 
wall in Fig. 11(a) is particularly striking. Figure 11 
(b) shows the difference between the initial state in 
Fig. 8(f) and the final relaxed state in Fig. 1 1(a). using 
a different color coding scheme. What is depicted 
here is the dot product between the magnetization 
directions of the initial and final states; each pixel 
represents a scalar value in the interval [ - 1. - [ ] 
corresponding to the dot product of the two unit 
vectors. In this coloring scheme * l corresponds to 
purple. - 1 to red. and the values in between are 
mapped continuously to the color wheel of Fig 1, 
Figure 1 1(b) clearly shows the nonuniform motion of 
the wall during zero-field relaxation. While some 
regions of the wall have not moved at all. others have 
travelled by as much as a few hundred angstroms. In 
general, the wall seems to have remained pinned on 
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several defects, while regions of the wall in between 
those defects have relaxed toward a state of mini- 
mum energy Once again, the observed behavior 
confirms that small patches of material with in-plane 
amsotropv are effective in capturing and stabilizing 
the small domains 

3. Concluding Remarks 

Hypothetical mechanisms of wall coercivity in 
thin films of amorphous rare earth-transition metal 
allovs were examined in this paper Lsing computer 
simulations, we found that regions as small as a few 
hundred angstroms :n diameter, with certain devia- 
tions of their structure magnetic attributes from the 
rest of the film, could trap domain walls and cause 
significant changes in the value of the coercive held. 
Values of the coercive field obtained by simulation 
are comparable to those observed in practice 
Whether or not these hypothetical sources exist in 
real materials is a question whose answer must await 
further progress in experimental nanomagnetics . 
Among the existing tools for observ ation of the mag- 
netic state in thin films. Lorentz Electron 
Microscopy 40 ’ and Magnetic Force Microscopy 411 
have the potential to clarify this situation in the near 
future. 
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